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Spin- and energy-dependent tunneling through a single molecule with intramolecular 

spatial resolution 
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We investigate the spin- and energy dependent tunneling through a single organic molecule (CoPc) 
adsorbed on a ferromagnetic Fe thin film, spatially resolved by low-temperature spin-polarized scan- 
ning tunneling microscopy. Interestingly, the metal ion as well as the organic ligand show a sig- 
nificant spin-dependence of tunneling current flow. State-of-the-art ah initio calculations including 
also van-der-Waals interactions reveal a strong hybridization of molecular orbitals and surface ?>d 
states. The molecule is anionic due to a transfer of one electron, resulting in a non-magnetic (S= 0) 
state. Nevertheless, tunneling through the molecule exhibits a pronounced spin-dependence due to 
spin-split molecule-surface hybrid states. 

PACS numbers: 68.37.Ef, 31.15.E-, 75.70.-i, 75.50.Xx 
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Molecular based systems are fascinating, yet promis- 
ing candidates for nanoscale spintronic devices and open 
viable routes toward quantum computing [1^. Previous 
experiments on the spin transport through break junc- 
tions [2] and spin valves [3 unveil exciting new frontiers of 
molecular magnetism. Much effort is dedicated to under- 
stand the properties of organic/magnetic interfaces. To 
this end, spatially averaging techniques show substantial 
spin-injection into, as well as long spin-coherent trans- 
port throughout films of metal-phthalocyanines (MPc) 
[Hl^. However, detailed and quantitative access to differ- 
ent constituents of a single molecule is desirable, though 
challenging. Scanning tunneling microscopy (STM) is 
well established as a probe of a local spin jHHTS] in an 
atomically well defined environment. 

lacovita et al. recently performed a spin-polarized 
STM (SP-STM) study of a CoPc in contact with a fer- 
romagnetic cobalt nano-island [TT. Stacking contrast, 
spin-dependent scattering, edge states, mesoscopic relax- 
ations as well as the adsorbate induced modification cre- 
ate a complex environment [15 toward understanding the 
influence of the substrate on molecular magnetism. Af- 
ter careful selection of electronically equivalent Co nano- 
islands a ferromagnetic exchange interaction between the 
molecular spin and the cobalt lead was successfully de- 
duced, both theoretically and experimentally. 

In this letter we demonstrate a significant spin- 
polarization for a CoPc molecule in contact with a fer- 
romagnetic Fe thin film due to molecule-substrate hy- 
bridization even though the molecule loses its net spin. 
As confirmed by SP-STM, an energy/site-dependent spin 
polarization from inversion to amplification is resolved on 
the sub- molecular scale. State-of-the-art density func- 
tional theory (DFT), which includes the decisive role of 



van-der-Waals (vdW) interactions, reveals both the mag- 
netic and electronic nature of the molecule coupled to 
the ferromagnetic substrate. Even though the net spin 
of the molecule is lost due to a transfer of one electron, 
spin-splitting is recovered through the local bonding of 
molecular orbitals with Fe 3(i bands. 

Simulations were carried out in the DFT [16 formal- 
ism with a plane wave implementation as provided by 
the VASP code |17j. Pseudopotentials used were gener- 
ated with the projector augmented wave method [18 by 
using the PEE generalized-gradient exchange-correlation 
energy functional [19^ (GGA). A slab consisting of two Fe 
and three W atomic layers, with a (5 x 7) in-plane surface 
unit cell modeled the molecule-surface system. The ki- 
netic energy cutoff of the plane waves was set to 500 eV 
while the Brillouin zone was sampled by the F point. 
Optimized molecule-surface geometries were obtained by 
relaxing all molecular degrees of freedom and those of the 
Fe overlayers by including long-range vdW interactions 
in a semi-empirical way [20l [21]. The threshold of the 
calculated forces was fixed to 0.001 eV/A. 

Experiments were performed in an ultra-high vacuum 
system with an STM operated at 6 K [22 . Chromium 
coated tungsten tips with out-of-plane spin sensitivity 
^23j were utilized for all measurements. 1.8 ML Fe was 
deposited on a W(llO) single crystal kept at 400K. After- 
ward, the magnetic sample was cooled down to measure- 
ment temperature and spin sensitivity was confirmed by 
imaging the alternating spin- up (0) and spin-down ((g)) 
domains of the 2^^ layer portion of the Fe film [23 . CoPc 
molecules were thermally sublimated from a home-built 
crucible onto the pre-cooled substrate to ensure spatially 
well separated molecules. 

All STM images were recorded in constant-current 




FIG. 1. (a) SP-STM image of CoPcs on magnetic iron do- 
mains (0,0). (b) Line profile as indicated in (a), (d)-(i) 
Higfi-resolution images of CoPcs in all three orientations (A, 
B, C) on oppositely magnetized domains, (c) Adsorption con- 
figuration of a CoPc in orientation A. (a) and (d)-(i) 0.05 V, 
200 pA; a) 240 A x 280 A; (d)-(i) 25 A x 25 A. The colorscale 
for (a) is non-linear to illustrate the domain structure, while 
the colorscale for (d)-(i) is linear. 



mode at the set current / and sample bias voltage U. The 
local spin-polarization of the atomic-scale tunnel junction 
can directly be determined from height variations As re- 
sulting from different magnetization directions of elec- 
tronically identical areas as introduced by Wiesendanger 
et al. ^: 



As<l A 
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The polarization of tip (Pt) and sample (Pg) are defined 
as the normalized difference of the energy integrated spin- 
up and spin-down states (Pt 
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^tfs) "-'-)• ^ ^^ ^^^ angle between the respective mag- 
netization directions, (j) the mean local tunneling barrier 
height, and a = 1 eV"^A" . 

Figure ll[a) shows a representative 2^^ ML area of Fe 



after deposition of CoPc molecules. The two dislocation 
lines indicate the crystallographic axes. Isolated CoPcs 
adsorb in three well-defined orientations (denoted as A, 
B, and C). The line profile (Fig. [lib)) along the indicated 
positions in (a) quantitatively illustrates the height vari- 
ation within the atomically fiat Fe film (As = 0.3 A) 
and the domain-dependent molecular appearance across 
a type A CoPc on an up- and down-domain. This vari- 
ation is due to the relative alignment of tip and sam- 
ple magnetization being either parallel or anti-parallel. 
While the center of the molecule (Co-site) has the same 
apparent height on both domains {As = A), we see 
a profound contrast at the ligand-site (As = —0.3 A). 
For an up-domain the ligand exhibits an increased ap- 
parent height relative to the surrounding Fe film, which 
is in contrast to the apparent height for a molecule on 
the down-domain. 

Thereby, orientation A [25 represents a highly sym- 
metric configuration (C2v-symmetry) with one molecular 
axis oriented exactly along the <001> direction. While 
there are differences in intramolecular contrast amongst 
the three molecular orientations (Fig. [ijd) to (i)), all 
three configurations show similar domain dependent con- 
trast variations. In the following, we focus on molecules 
in configuration A which are not in the vicinity of disloca- 
tion lines or domain walls. As the molecular appearance 
controllably switches under an external magnetic field 
applied [26 , the magnetic contrast can unambiguously 
be attributed to the spin-polarized local density of states 
(SP-LDOS) of the molecule-surface system. 

First-principles calculations within DFT (GGA) clarify 
the origin of the observed spin-dependence. Calculations 
for the A geometry show the preferred top adsorption, 
i.e. the Co atom on top of a surface Fe atom. Previ- 
ous results for a Co island revealed an adsorption at the 
bridge site [14 , while the hollow site was deduced for a 
Au(lll) substrate [12 . Furthermore, calculations cover 
both spin-orbit coupling (SOC) as well as the consider- 
able role of vdW interactions. Results without SOC are 
discussed in terms of molecular orbitals (MOs) to clearly 
illustrate the effect of vdW interactions alone. 

Figure [2] depicts the change in electronic structure of 
a CoPc with increasing molecule-substrate interaction. 
Figure [2|a) shows the SP-LDOS for the free molecule. In 
agreement with previous work [9l [121 113 [21] the origin of 
spin-splitting is an unpaired electron in a MO with d^^ 
contribution situated at the Co-site and the total molec- 
ular spin S is 1/2. Figure [2|b) presents the effect of the 
surface, when vdW forces are neglected during the relax- 
ation process. The molecule adsorbs 3.1 A above the sur- 
face remaining fiat. The SP-LDOS shows a hybridization 
of MOs and substrate 3d states. According to the spa- 
tial extension of the MOs perpendicular to the molecular 
plane (i.e., with tt character), the hybridization of MOs 
containing the Co dz2 atomic orbital is the strongest, fol- 
lowed by those including d-j^ and Pz atomic contributions. 
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FIG. 2. The geometry and electronic structure for a free CoPc 
(a) and adsorbed on an Fe surface (b) without and (c) with 
vdW forces included during the relaxation, (d) Molecular 
deformation due to vdW forces. The Ce rings and two of the 
outer N are twisted by up to 0.3 A toward the surface and 
away from the plane defined by the Co ion and four inner N. 
All H are pointing away from the surface. 



The da- and a MOs are only slightly broadened, com- 
pared to the free molecule case, as they are localized in 
the molecular plane. The spin-splitting of the molecule- 
surface hybrid states is reduced due to a transfer of an 
electron from the substrate to the CoPc. As a result the 
formerly unoccupied d^^ type MO becomes occupied and 
the total molecular spin is quenched (0 fib)- 

The role of vdW forces is crucial as it brings the 
molecule 0.5 A closer to the surface (Fig. M^b)) and dis- 
torts the molecular geometry (Fig. [2|d)). This new ad- 
sorption geometry has a drastically different electronic 
structure due to the overall changes of hybridization be- 
tween molecule and substrate. The (SP-LDOS Fig. [2|c)) 
shows not only MOs with a tt character strongly hy- 
bridize with spin-polarized Fe 3d states of the same sym- 
metry to form broad spin-split bands but also a type 
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FIG. 3. Experimental and simulated SP-STM images at 
U = 0.05 V and 22 A x 22 A for both spin directions, (a) Av- 
eraged experimental images, (b)-(d) Isocharge surfaces; (b) 
conventional DFT (GGA) (without vdW), (c) DFT+vdW, 
(d) DFT+vdW and including SOC. 



MOs are significantly affected by the interaction with 
the surface. Again, a transfer of one electron from the 
surface to the molecule annihilates the molecular spin, 
but spin-splitting is recovered due to the local molecule- 
surface bonding at different parts of the molecule. The 
newly formed molecule-surface hybrid states have, within 
a given energy interval [E^E^], an unbalanced, locally 
varying electronic charge in the up and down channels 
which is mapped in SP-STM {E = Ef,E' = Ef^ eU). 

For a direct comparison of the first principles calcula- 
tions with constant-current SP-STM images, isocharge 
surfaces above the CoPc are extracted from the spa- 
tial variation of the energy integrated SP-LDOS. This 
approach mimics the experimental situation of a local 
and spin-sensitive tip probing the charge density above 
the surface at constant current and thereby, accounts 
for the variation of decay lengths of the different states 
into vacuum. Figure [3] gathers calculated spin dependent 
isocharge surfaces and experimental SP-STM images for 
U = 0.05 V. The isocharge surfaces reproduce all fea- 
tures of the experimental SP-STM images, only if vdW 
forces forces are taken into account. Corrections due to 
the inclusion of spin-orbit coupling effects are minor. 

Figure [4] illustrates quantitatively the local spin polar- 
ization. To compensate for the tip polarization, we nor- 
malized the experimental polarization above the molecule 
(^CoPc) i"6lative to the spin polarization of the free sur- 
face (Ppe^), 



pexp 
^CoPc 
pexp 



PtPi 



CoPc 



cosO 



PtPpeCOS^ 



TDsim 
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The normalized polarization (P) is therefore, to a first 
approximation, independent of the tip polarization and 
can directly be compared to our simulated polarization 

/ psim\ 

For comparison both experimental and theoretical spa- 
tial maps of the polarization at three characteristic ener- 
gies are given in Fig. [4j We see an excellent agreement 
of experimental findings and DFT calculations, when 
the vdW relaxed adsorption geometry is used. Remark- 
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FIG. 4. Local spin polarization at (a)-(c) three representative 
energies. Raw- and avg (over multiple images and according 
to symmetry [26 ) experimental data are compared with DFT 
simulations inch vdW and inch Spin-Orbit-Coupling. Line 
profiles follow high-symmetry directions as indicated in the 
sphere model inset in (a). Simulated data including SOC is 
only given in (b) as the SOC corrections for (a) and (c) are 
negligible. Insets: 22 A x 22 A. 



ably, not only reduction (1 > P > 0) but also inversion 
(P < 0) and amplification (P > 1) of spin polarization 
can equally be observed. The P above the Ce-ring shows 
an inversion for all energies, while the spin polarization 
near the Co-site ranges from inversion (at -0.2 V) to am- 
plification (at 0.75 V). 

In conclusion we investigated single CoPc molecules 
adsorbed on a ferromagnetic Fe surface and locally ob- 
served the energy-dependent spin transport with in- 
tramolecular spatial resolution. The tunneling current 
above the entire molecule - metal ion, as well as the or- 
ganic ligand - shows a high, locally varying spin polariza- 
tion ranging from inversion up to amplification with re- 
spect to the ferromagnetic Fe film. Optimized molecule- 
surface geometries for the first principles calculations 
were obtained by including long-range van-der-Waals in- 
teractions during the relaxation, leading to an excellent 
agreement of theoretical and experimental data. There- 
fore, the observed spin polarization is identified as a 
unique property the molecule-substrate hybrid states, by 
our combined DFT and SP-STM approach. 
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